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fc-core percolation is a percolation model which gives a notion of network functionality and has 
many applications in network science. In analysing the resilience of a network under random damage, 
an extension of this model is introduced, allowing different vertices to have their own degree of 
resilience. This extension is named heterogeneous fc-core percolation and it is characterized by several 
interesting critical phenomena. Here we analytically investigate binary mixtures in a wide class of 
configuration model networks and categorize the different critical phenomena which may occur. Wc 
observe the presence of critical and tricritical points and give a general criterion for the occurrence 
of a tricritical point. The calculated critical exponents show cases in which the model belongs to the 
same universality class of facilitated spin models studied in the context of the glass transition. 



I. INTRODUCTION 

A significant issue in many networked infrastructures is 
the threat that local damage can represent to the whole 
system [TJH]. Moreover, there is an important distinc- 
tion in the way a network collapses: it can happen in 
a smooth, predictable way, or with extreme abruptness. 
The capability of designing networks where a collapse can 
be predicted in advance is expected to be increasingly 
important [5]. 

In recent years, there has been a wide renewal of inter- 
est in percolation models, regarding the occurrence of a 
first order transition, instead of the classical continuous 
percolation transition [3]. However, it has been also shown 
[5 that discontinuous transitions are not so unusual in 
percolation models and can be dated to severals years 
ago [HIE]- Moreover, several examples have been recently 
proposed in the context of explosive percolation [8HTT]. 
interdependent networks [12] and hierarchical lattices [13] . 

fc-core percolation is an extension of percolation provid- 
ing in a simple model a wide range of critical phenomena: 
a giant fc-core cluster may collapse either continuously 
or discontinuously as a function of random damage [TH - 
[15] . Moreover, the analysis of the fc-core architecture of 
a network has developed into applications in different 
areas of science including protein interaction networks 
[17] . jamming [13], neural networks [15], granular gases 
[20] . evolution |21j . social sciences [32] and the metal- 
insulator transition [23]. Finally, a set of spin model 
approaches to the glass transition shares close similarities 
with fc-core percolation [24]. 

The analytical formalism of fc-core percolation has been 
recently extended to include a local notion of robustness: 
some nodes can be more resilient than others and require 
a smaller number of neighbors to remain active. This ex- 
tension has been named heterogeneous fc-core (HKC) per- 
colation and has been studied in locally tree-like networks 
|25j . finding a number of interesting critical phenomena 
including a tricritical point (TCP) [26] . 

In this paper, we extend the analysis of [35] to present 



an exhaustive description of binary mixtures in hetero- 
geneous fc-core percolation and show that heterogeneous 
fc-core percolation models are characterized by a wealth 
of critical behaviors, involving both first and second or- 
der transitions. In Section 2 we define heterogeneous 
fc-core percolation and explain the mathematical formal- 
ism in locally tree-like networks. In Section 3 we focus on 
some illustrative examples of binary mixtures of vertex 
types, examining the different critical phenomena and 
calculating relevant critical exponents. In Section 4 we 
give a general argument of the occurrence of a tricritical 
point in a binary mixtures of vertex types in heteroge- 
neous fc-core percolation. Finally, Section 5 states our 
conclusions. 



II. THE MODEL 

Given a simple graph, a fc-core is defined as the largest 
subgraph where every vertex has at least fc neighbors in 
the subgraph itself. A common way of determining the fc- 
core consists in removing recursively all the vertices (and 
adjacent edges) with less than fc neighbors. We consider a 
framework where the nodes (and adjacent edges) of a given 
network are randomly removed with probability 1 — p and 
we ask how the fraction of vertices in the fc-core varies as p 
is decreased. An analytical formalism has been developed 
to study this problem on the configuration model [16] . The 
configuration model is defined as the maximally random 
network with a given degree distribution P(k). It has the 
important property that the number of loops vanishes 
as the size N — > oo, which guarantees that if a fc-core 
exists, it must be infinite, at least if fc ^ 2 [151 H3- This 
formalism is based on the definition of (fc — l)-ary subtree. 
Given the end of an edge, a (fc — l)-ary subtree is defined 
as the tree where, as we traverse it, each vertex has 
at least fc — 1 outgoing edges, apart from the one we 
came in. In configuration model networks, then, the fc- 
core coincides with the (fc — l)-ary subtree. The formalism 
by Dorogovtsev et al. solves the problem in any locally 
tree-like graphs and determines the order of the transition 
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with which the giant fc-core collapses [16] . 

The heterogeneous fc-core is an extension of the fc- 
core where the minimum threshold fcj is a local parameter, 
i. e. it depends on the vertex i. The (fcj — l)-ary subtree, 
then, is the tree in which, as we traverse it, each encoun- 
tered vertex has at least fcj — 1 child edges. We define Z 
as the probability that a randomly chosen vertex is the 
root of a (fcj — l)-ary subtree and M as the probability 
that a randomly chosen vertex is in the HKC. We can 
then write M ab for a mixture of two types of vertices 
k a , k b where k a and kb vertices occur with probability r 
and 1 — r, respectively: 

M ab (p) = M a (p) + M b (p) (1) 

oo k a 
q—k a l—q 

oo kb 
q=k b l=q 

where M a i b \[p) is the fraction of nodes of type a(b) in the 
heterogeneous fc-core, respectively, P(q) is the degree dis- 
tribution of the network and we have used the convenient 
auxiliary function: 

*i, 9 (*, z) = (fj (i - zf- x ( 1 )x m (z - xy- m . 

^ ' m=l ^ TO ' 

It is shown in that Z satisfies the self-consistent 
equation: 

q=k a W l=q-l 

oo p( \ ^'fJ — 1 

E ( 4 ) 

q=k b W l=q-l 

In general, it is not granted that the heterogeneous 
fc-core is uniquely made up of a giant cluster. As it is 
possible that some finite clusters of the heterogeneous 
fc-core are present, Baxter et al. develop the equation for 
X, the probability that an arbitrarily chosen edge leads 
to a vertex which is the root of an infinite {ki — l)-ary 
subtree. In the case of a binary mixture, the equation 
reads 

q=k a W l=q-l 

oo p( \ fcfa — 1 

+p^-r) z Z q rr £ ( 5 ) 

q=k b W l=q-l 

Therefore, the probability S ab that a randomly chosen 
vertex belongs to the giant HKC is 

S ab = £ 5 «(0 + X) S 6 (Z) ( 6 ) 

;>fc Q z>fct, 



with 

5 o (0 = pt^;p(«)$ I(9 (a: ) z), (7) 

^((1=1,(1-^^(1,2). (8) 

?>i 

An important concept in fc-core percolation is the 
corona. The corona is defined as the subset of the HKC 
where every vertex i has exactly ki nearest neighbours in 
the HKC. By definition, then, a corona cluster is char- 
acterized by the property that if only one of its vertices 
is removed, the whole cluster collapses. It can be shown 
that the corona clusters are finite everywhere except at 
the phase transition, where the mean cluster size diverges 
[231 12"5] . Using equation ([6| , then, the corona C ab of a 
binary mixture is given by the following formula 

C ab = S a (k a ) + S b (k b ) = 

= P rY P{q)$k a , q ( X ' Z ) + 
q>k a 

+p{l-r)Y J P (q)^q{X,Z). (9) 

q>k b 

The distinction between heterogeneous fc-core clusters 
and the giant heterogeneous fc-core has to be made when- 
ever the two do not coincide. That is the case when we 
consider binary mixtures of the type k a = 1, k b > 1, as 
for k = 1 finite 1-core clusters are possible. However, for 
mixtures such as 2 < k a < k b , no finite heterogeneous 
fc-core are possible, therefore: X = Z and M ab = S ab . 

As in the homogeneous case, also in heterogeneous fc- 
core percolation a fc-core architecture of the network can 
be defined. Therefore we define a heterogeneous sub-core 
as a subset of the HKC where a higher threshold ki is im- 
posed on some or all vertex types. More specifically, given 
a HKC of type k a , k b , the strength Mj~ of a heterogeneous 
sub-core is given by 

M k {p) = M% + M\ = 

= r P Y P(q)$h a Az,z) + 

q>ha. 

+(l-r)pJ2 P (<l)^h b , q (Z,Z), (10) 

q>h b 

where h a = max(k. k a ) and h b = max(k, k b ). An analo- 
gous expression holds for 5fc(p). Therefore, when consid- 
ering for example a binary HKC with k a < k b , sub-cores 
with k < k a coincide with the HKC, whereas sub-cores 
with k > k a are proper subsets of the HKC. 

III. ILLUSTRATIVE EXAMPLES 

We consider now a few examples of binary mixtures 
on Erdos-Renyi graphs with the most interesting critical 
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phenomena. Plugging in the Poissonian degree distribu- 
tion P(q) = z\ exp(— Z\)/q\ (where z\ is the mean degree) 
into equations Q and ([5]), the sums can be calculated 
analytically and several critical phenomena can be found. 
First, we review two cases already studied in the liter- 
ature: the case k = (k a ,ki,) = (1,3), examined in |25j 
on the Bethe lattice, and the case k = (2,3) [25]. Then, 
we focus on the case k = (2,4), and show that its phase 
diagram is qualitatively identical to the ones of the type 
k a = 1, fcb > 3, and represents a quite general behavior, 
with applications to the theory of glass transitions. 



Using ([9]), the strength of the HKC corona is given by the 
following expression: 



pe' ZlZ [rziX- 
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(1 - r)zf (3Z 2 X - 3ZX 2 + X 3 



(20) 



In particular, the fraction of vertices of each type in the 
HKC are given by the following expressions: 



Si = pr(l 



,-ziX 



(21) 



The case k„ 



1, k b 



In this case we have to take into consideration that, 
due to the presence of vertices of type 1, there are finite 
A:-cores. The two equations Q and ^ can be re-written 
as 



where 



hs(Z) 



pfu(Z) = 1 
ph 13 (X, Z) = l, 

1- (l-r)e- ZlZ (l + z-lZ) 



h 13 (X, Z) = 1 - (1 - r)z t < 



— z\Z 



(11) 

(12) 



(13) 



(14) 



It can be shown that the locus defined by f[ 3 (Z) = 
corresponds to a line of first order transitions which ends 
into a critical point defined by f[ 3 (Z) = f" 3 (Z) = [29] . 
A simple calculation yields: 



1 

r c = 1 e 

3 



3 

Z\ 



Z r 



1 

Z\ 



(15) 
(16) 
(17) 



X c has no analytic expression, being the non-trivial solu- 
tion of the equation 1 — exp(— ziX c ) = 2/3ziX c , which 
is ziX c ~ 0.874. 

Using ([6| , we can calculate the strength of the giant 
HKC as 

S 13 =p[l~e-^ x ] + 

-p(l - r) Zl Xe-^ z (l + z x Z - \zix\ , (18) 



and we can also derive the expression of the 2-sub-core: 
S 2 = pr{l - e~ ZlX - v^-ziZ-s 
+p(l -r)[l-< 



z\Xe 

-ZxX 



-Z\Xe 



1 



l + Zl Z-- Zl X 



(19) 



S 3 =p(l-r) 



1 - e 



-ziX 



z\Xe 



1 + ZiZ- \z x X 



(22) 

The critical percolating strengths can be calculated 
with arbitrary precision. The approximate values are 
given in Table I. 



r c ziSfz Z1C13 Z1S1 Z103 
0.0939 0.3821 0.2700 0.1642 0.2179 

TABLE I. Critical fractions of some relevant quantities in the 
case k = (1, 3). 

Figure [T] displays the phase diagram of the mixture 
k = (1, 3) for a Erdos-Renyi graph. The two phase region 
occurs at a relatively low value of r c , meaning that many 
3-nodes are necessary to drive the system towards the first 
order (hybrid) transition. The relative composition of the 
giant HKC, instead, presents a much higher fraction of 
1-nodes, due to the high fragility of 3-nodes with respect 
to 1-nodes. The phase diagram also shows that the corona 
strength is much smaller than S± 3 on the right side of the 
coexistence region, whereas is much closer to Si 3 on the 
left side, showing that the corona clusters dominate the 
HKC in the 1-rich phase. 



B. The case k a = 2, kt = 3 

As shown in our recent paper, this case has the property 
X — Z and is characterized by a tricritical point [26] . We 
review the peculiarities of this case. The quantity Z at a 
given damage fraction (1—p) can be calculated by solving 
the equation pf2 3 (Z) = 1, where 



l-e- ZlZ [l + (l-r) Zl Z] 



(23) 



From 0, the strength of the heterogeneous fc-core is 
given by 



M 23 (p)=p{l 



- Zl z 



1 + z x Z + -(1 - r)z\Z 2 



(24) 

It is easy to see that f2z(Z) has a maximum at a finite Z 
for r < 1/2, which continuously moves to Z — at exactly 
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FIG. 1. Phase diagram of the case k = (1,3) for the Erdos- 
Renyi graph (z\ = 10). The main panel illustrates the critical 
values of the giant HKC strength. The blue lines report the 
second order (solid) and the first order (dashed) transitions. 
The critical (black) and the critical end point (blue) are indi- 
cated by a dot. The fractions of nodes of type 1 (<Si) and 3 
(S3) in the giant HKC are also reported. The corona size C13 
is represented by a magenta dashed line. The inset reports the 
most relevant detail of the phase diagram in the space (r,p c ). 



for r < 1/2, which smoothly shrinks at the TCP. The 
fractions M 2 and M 3 of nodes of each type inside the 
HKC are also plotted. It is interesting to note that M 2 is 
substantially higher than M3 in the proximity of the TCP. 
It appears that the vanishing of the discontinuity is finely 
tuned by the vanishing fraction of type 3 nodes in the 
HKC. At the same time, of course, the fraction of type 2 
nodes becomes increasingly important, determining the 
change of order of the transition at the TCP. 
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r = 1/2. This implies that there is a line of first order 
transitions for r < 1/2, and a second order percolating 
transition for r > 1/2. At r = 1/2, the two lines match 
exactly at a tricritical point (TCP). The critical exponent 
/3 at the transition M^^ip) ~ (p~Pc)^, can be calculated 
analytically: 



(25) 



The value of the exponent /3 at r < 1/2 agrees with the 
typical hybrid transition phenomenology |15[ 116] , whereas 
for r > 1/2 we recover the exponent of classical percola- 
tion without dangling ends [33]. The value of j3 at the 
TCP does not change when it is calculated along a line 
at fixed p c . 

Another exponent which can be calculated analytically 
is the one which governs the vanishing of the coexistence 
region as r — > l/2~. We indicate it as M*(r) ~ (| — r)' 3 ", 
with /3 U = 2. Finally, the rotation defining the critical 
fields is 




cos 9 sin 9 
— sin 9 cos 




(26) 



with tan6> = A/z\. Close to the tricritical point, the 

1/2 

critical line has a behavior m ~ , with a crossover 
exponent tpt = 2. 

Fig. [2] shows the phase diagram of this case. The 
system undergoes a first order (hybrid) phase transition 



FIG. 2. Phase diagram of the k = (2,3) mixture, showing 
the total mass of the percolating heterogeneous fc-core cluster 
at different compositions r, for Erdos-Renyi networks with 
z\ — 10. The TCP at r = 1/2 separates a line of first order 
transitions (dashed) from the second order line (solid). The 
masses of the fractions of nodes of types 2 (green) and 3 (red) 
in the giant HKC are also shown. The inset shows the phase 
diagram in the (r, p) space. 



C. The case k a — 2, k b — 4 

From equation Q, the case k a ~2. fc& = 4 is solved by 
Phi{ z ) = 1> where 

f u{ z) = i- e -' lZ {i + (i-0(^ + l^ 2 )} . (27) 



Differently from the case k a — 2, kb — 3, here for some 
values of r the function / 2 4 has a local maximum which 
does not approach the Z = point. More specifically, it 
is easy to show that fiA{Z) has three types of behavior: 
monotonously decreasing, a local maximum at Z > 0, and 
a global maximum at Z > (Fig. [3]). A global maximum 
at Z = corresponds to a second order de-percolating 
transition, as the solution of equation f(Z)p = 1 smoothly 
vanishes at some p c . A local maximum at Z > 0, instead, 
determines a first order transition line between two stable 
giant HKCs. 

Using equation ([lj, the strength of the HKC can be 
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(a)r = 



(b)r = 0.15 




(d)r = 0.25 



(e)r : 



(f)r = 0.35 



FIG. 3. Different scenarios of the solutions of equation 
pf2n{Z) = 1 in the case k a = 2, kt = 4 for different val- 
ues of r (zi = 10). The red horizontal line marks the value 
of 1 /p c where the first order transition occurs. In the panels 
(a) and (b) there is a first order transition corresponding to 
the global maximum of f2±{Z); in (c) the maximum is at the 
same height as ji^iZ = 0) (critical end point); in (d) there 
is both a first and a second order transition; in (e) the first 
order transition disappears into a critical point; in (f) only 
the classical percolation transition at Z — remains. 



re-written as 

M 24 (p) =J5 jl-e- 2lZ 
1 



l + Zl Z+-(l-r)zlZ 2 - 



+ ^(l-r)z 3 Z 3 



(28) 



From the definition (|10|, the strength of the heterogeneous 

(29) 



sub-core as a subset of the HKC is given by 

M k (p)=Ml+M b k , 
where 

ha-l 



M a k =pr 



M b k =p(l-r) 



{ Zx zy 



(30) 



n=0 



n 



and h a (b) = max(k, k a n,)). For k < 2, the sub-core coin- 
cides with the HKC and for instance we have Mi = M24. 
For k > 4, the sub-core is a proper subset of the HKC, 
but the vertex distinction is not longer relevant, as the 
threshold k applies in the same way to both vertex types. 
Therefore, the sub-core coincides with the usual homoge- 
neous fc-core. The most interesting case is when k — 3, 
as this threshold restricts the number of acceptable type 
4 vertices, but it has no effect on the threshold of type 2 
vertices. We indicate the corresponding strength as M23, 
which is given by the following formula from (29): 



M 23 (p) =p\l 



-Z\ Z 



1 



2 rj2 



1 + z\Z + -zfZ 



r)zfZ 3 



(32) 



From (27), we calculate the phase diagram at different 
compositions r (Fig. |4|. At high r, the phase diagram is 
characterized by a critical line of de-percolating transi- 
tions. This line meets a first order line at a point which 
is usually called a critical end point in condensed matter 
physics. The first order line corresponds to the type of 
hybrid transition observed in fc-core percolation for fc > 3. 
Differently from the case k a = 2, fc& = 3, here the first or- 
der line presents a critical point at the end of a two-phase 
coexistence between a low and a high density phase. 
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FIG. 4. Phase diagram of the case k = (2, 4) for the Erdos- 
Renyi graph (zi — 10) in the space (r,M), main panel, and 
(r,p), inset. Symbols are as in Fig. [l] 



The position of the critical end point is determined by 
imposing fz^Z) = and f2i(Z) — rz±, which yield the 



soluti 



0.2239964831566. Now we consider the 



strength M| 4 (r) of the HKC along the low density border 
of the coexistence region in approaching the critical end 
point. We define the critical exponent f3 cep of the critical 
(31^ end point from the manner in which such low density 
border vanishes: M| 4 (r) ~ (r — r cep ) l3 " p . It emerges that 



G 



as r — > r2 



we have: 



(33) 



which implies that f3 cep — 2. The strength of the 3-sub- 
core in the HKC can be calculated in a similar way 



M*{r)~{r-r cep f. 



(34) 



The higher exponent of M| with respect to implies 
that the critical 2-nodes are essential for observing the 
critical end point. The k = 3 threshold, in fact, recursively 
eliminates all the type 2 nodes with exactly two neighbors, 
triggering a cascade which essentially makes the sub-core 
become a negligible fraction of the HKC. 

The critical point can be calculated by requiring that 
the maximum (/^(Z) = 0) coincides with the second 
change in convexity (f^iZ) = 0). Hence we get an 
equation for the critical composition r c : 



exp 



(l + \Iy^t) = 5 _ 7fc + V(l-r c )(l-3r c ), 

(35) 

with numerical solution r c w 0.2962590188. Using the 
appropriate equations we can calculate with arbitrary 
precision the critical damage p c and the values of the 
k-cove strengths: they are summarized in Table II. 



Z!p c 2iAf 2 C 4 2lC| 4 ZlMf Z\M% Z 1 M§ 



0.2963 4.1131 0.6510 0.4079 0.4963 0.1547 0.3569 

TABLE II. Critical fractions of some relevant quantities in the 
case k = (2, 4). 

The critical exponent j3 defined by M24 ~ (p — p c )^ 
can be calculated analytically for each region of the phase 
diagram: 



r > r, 



2 critical line 
1/2 hybrid transition r <r c 
1/3 hybrid transition r = r c 



(36) 



In the region r cep < r < r c , the first order transition sep- 
arates two percolating phases and therefore it is possible 
to calculate the exponent /?' on the left hand side of the 
first order line (Fig. El): 



1 

1/3 r = 



(37) 



The exponent /3' = 1 away from the critical point is not 
singular, whereas from the other side of the first order 
transition it is j3 — 1/2. This difference is due to the 
presence of the hybrid transition, which is asymmetric. 
However, the two exponents coincide at the critical point 
as generally expected. 

It is interesting to note that the exponents at the hybrid 
transition in ( |36[ ) exactly match the ones found in facili- 
tated spin models reproducing mode-coupling theory sin- 
gularities [3D] . In particular, exponent 1 /2 corresponds to 



the Ai singularity associated with a discontinuous liquid- 
glass transition, whereas the exponent 1 /3 corresponds to 
the A 3 singularity at the endpoint of the discontinuous 
glass-glass transition of the F13 schematic model [3TH 131j . 

The difference AM24 of the two coexisting HKC 
strengths as r — > r~ vanishes according to the follow- 
ing law: 



AM 24 ~ (r c - r) 



1/2 



(38) 



Therefore, the associated critical exponent is j3 u = 1/2. 
The expansion of AM3 yields the same critical exponent 



AMs-^-r) 1 / 2 . 



(39) 



The phase diagram, then, is characterized by the same 
topology as in the k a = 1, k{, = 3 case: a first order 
(hybrid) transition line which ends in a critical point and 
a critical line which encounters the first order line at a 
critical end point. Here the critical point occurs at a 
higher fraction of low-fc nodes, and at a larger fraction of 
low k nodes in the composition of the HKC: the opposite 
of the case k = (1,3). This is presumably due to the 
higher fragility of 4-nodes with respect to 3-nodes. 



D. The case k a = 3, kt — 8 

Binary mixtures involving thresholds higher than 2 
cannot present continuous transitions. However, the mix- 
ing of different types of nodes may still result in a crit- 
ical point. Fig. [5] shows the phase diagram of the case 
k = (3, 8) as an example of this behavior. The high dif- 
ference in resilience between 3-nodes and 8-nodes results 
in an intermediate region of the phase diagram where an 
8-rich phase collapses into a 3-rich phase. The two lines 
of first order transitions do not match smoothly, and a 
region with a stable 3-rich percolating phase can be easily 
seen. It can be shown that the critical point is in the 
same universality class as the one in the case k = (2,4), 
with the same critical exponents along the first order line. 



E. Summary 

Fig. [6] summarizes the phase diagram of a few relevant 
binary mixtures. The general behavior of mixtures of type 
k = (1, k) or k = (2, k) is characterized by a critical line 
and a line of first order transitions which ends in a critical 
point. The case k = (2, 3), however, is quite peculiar as 
the two lines match at a tricritical point. This occurrence 
is not obvious and will be investigated in the next section. 
Both the tricritical and the critical point observed in 
most mixtures belong to the same universality class of 
facilitated spin models reproducing mode coupling theory 
singularities of type F 12 and -F13, respectively [301132] , 

Fig. [f| also reports the case k = (3, 4) as an example 
where there cannot be continuous transitions, as both 
values of k are larger than 2. The mixture k = (3, 4) is 
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FIG. 5. Phase diagram of the case k = (3, 8) for the Erdos- 
Renyi graph (zi = 30). Symbols are as in Fig. [I] Differently 
from the previous cases, here there are two lines of first order 
transitions and no classical continuous percolation scenario. 



IV. A CRITERION FOR A TRICRITICAL 
POINT 



We have already noted that the occurrence of a TCP is 
quite peculiar in this type of models, but quite resilient to 
different network topologies 26J. In this section we give 
some insight into the origin of a TCP and show a criterion 
to establish the occurrence of it in a given binary mixture 
of vertices. As a general rule, given a binary mixture 
k = (k a , kb) where k a - and fc^-cores are characterized by a 
second and a first order transition, respectively, the TCP 
occurs whenever the critical point associated with the 
discontinuous transition of the fc^-rich mixture coincides 
with a vanishing HKC strength. The critical percolation 
line, by definition, is characterized by a vanishing HKC 
and must match the first order line to give origin to a 
TCP. 

In the following we assume that the degree distribution 
P{q) satisfies 



only characterized by a line of first order transitions. As 
we have seen in the case k = (3, 8), though, a critical point 
can still arise when the two values of k are quite far from 
each other so that the high-fc phase can collapse without 
undermining the stability of a \ow-k phase (Fig. [5]). 




FIG. 6. Comparison of the phase diagrams of the cases 
k = (1,3), k = (2,3), k = (2,4) and k = (3,4) for Erdos- 
Renyi graphs with z\ — 10. As in the previous phase dia- 
grams, continuous lines represent continuous phase transitions, 
whereas dashed lines represent first order phase transitions. 
The red dot represents a TCP; the other dots are critical 
points. 



0< 1)( 9 -2)P( 9 ) <+oo. (40) 

9=3 



This condition is quite general, as it has already been 
shown that scale free networks with P(q) ~ q~ J (for large 
q) and 1 < 7 < 3 are not characterized by first order 
transitions and therefore their phase diagrams cannot 
involve neither critical or tricritical points |16| . Therefore, 
the only interesting case where the following analysis does 
not apply involves scale free networks with 3 < 7 < 4. 

Due to the different properties of binary mixtures, we 
differentiate two cases according to the presence of finite 
HKC clusters or not. 



A. The case k a = 2, kb = k > 3 

In the case k = (2, k), equation (Q can be re- written 
as pf2k(Z) = 1, where the relevant function is 

&m= r £^g( 5 7 1 ) z -'(i-^-'-' + 

q=2 W 1 = 1 ^ ' 

Q — fa / — — X 



In this case X = Z, so we only need to study the behaviour 
of the equation at Z — > + . The expansion reads 



q=2 



q(q-l)P(q) 



+ OO 

J r E 



(q) 

q(q-l)(q-2)P(q) 



2 £ <g) 



0{Z 2 



-(i->£^(r>'- 2 +°<^'>- 

(42) 



q—k 



A necessary condition for a critical point which contin- 
uously approaches the line Z = for r — > r c is that 



/2 fc (0) = at r = r c . From equation (42), we have two 
possibilities. If k > 4, 



/»(0) = -gr 



+ 00 

E 

q =2 



g(g-l)(g-2)P(g) 



(43) 



which is always negative and vanishes only for r = 0. If 
fe = 3, instead, the linear term proportional to r mixes 
with the linear term from the type 3 nodes and we get 



/2 3 (0) 



E 

9=3 



g(g-l)(g-2)P(g) 

(?) 



(44) 



which vanishes at r = 1/2. 



From equation (42) it transpires that Zf 2 k{Z) can be 



written as an expansion in terms of core shells for each 
node type. The condition of the critical point f^{Z) = 
implies that the dominant term for Z — > + is the lin- 
ear one. The term iii^ 1 )^-^) p ^i) £2 re p resen ts the 
probability that, following an edge leading to a node of 
type 2, there are two outgoing edges connected to the HKC 
(i. e. there must be exactly one extra neighbor beyond the 
minimum possible). The term (1 — r) ^2 q (feli)^ fe l 
represents the probability that, following an edge leading 
to a node of type k, there are k — 1 outgoing edges con- 
nected to the HKC. In other words, this is the probability 
that the fc-node found at the end of the edge is connected 
to the /c-corona. The condition for a TCP, therefore, is 
that the 3-shell term of the type 2 nodes has the same 
order of magnitude as the fc-corona of the k type nodes. 
This can only occur when k = 3, where we have a TCP 
at r = 1/2 (see also Eq. [44 ) . 



B. 



The case k a = 1, = k > 3 



As noticed in Section 3, finite fc-cores exist when k = 
1 [26 . Therefore, we have X < Z, and the functions 
associated with the two relevant equations Q and ^ 



+oo 



q=k W 



9-1 

E 

l=k-l 



q l ^z l -\\-zy- l -\ 



(45) 



h lk (x,z) = E ( q V" 1 " 'C 1 - *) a_1 - m + 

q=0 ro=l ^ ' 



(q) 

q=k N ' J=fc-1 



(46) 



As in the case k = (2,fc), it is meaningful to expand 
h(X,Z) for X -> 0+: 



h lk (X, Z) = a {Z) + ai {Z)X + 0(X 2 ), (47) 



where 



rt 



q=2 
9-1 



q—k 



(q) 



J2 (* ^yZ^Q.-Z)*- 1 - 1 ; (48) 



ai{Z) 



-, +OO 

rE 



g(g-l)(g-2)P(g) 

(q) 



2 ^ (g) ^ \ I 

q=k w i=fe-l 

a(i-i)z l - 2 (i- zy- 1 - 1 



(49) 



It is important to remark that a\{Z) < for every Z, 
i. e. /i(0, Z) is always a local maximum of h as a function 
of X. 

If a TCP exists, there must exist a value of r for which 
a critical point (the endpoint of a line of first order tran- 
sitions) occurs exactly at X = 0, as this corresponds to 
the critical line. As lim^^o fik(Z) = +oo for every r > 
(Eq. 45), the equation pfxh(Z) = 1 never has a solution 
Z = and fi k is monotonic in the neighborhood of Z = 0. 
Therefore, a line of first order transitions corresponds 
to a local maximum of f\ k , where equation pf\ k (Z) = 1 
has two solutions. The discontinuity in the solution Z(p) 
provokes a discontinuity in the solution X{p) of the sec- 
ond equation. This in turn causes a discontinuity in the 
strength S of the giant HKC. Therefore, if a critical point 
occurs at Z = Z c , it must be f' lk (Z c ) = f" k (Z c ) = 0. 
(The local maximum of / cannot occur at Z = 0, as we 
have seen.) In order to have a TCP, the above defined 
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critical point must occur at X c — 0. This implies that 
the fraction p c of undamaged nodes at criticality must 
satisfy the inequality 



h lk (Q,Z c ) < 



1 

Pc 



Using the first equation, this is equivalent to 
h lk (0,Z c )<f(Z c ), 



(50) 



(51) 



where f[ k (Z c ) = f" k (Z c ) — and, in case of multiple 
solutions, the value of Z c where f\k{Z c ) is highest must 
be considered. 



The condition (51) is only a necessary condition for a 



TCP, because in general it may be possible that h has a 
maximum higher than /iifc(0, Z). However, we now show 
that this condition is violated for all value sof k. Let us 
consider the two conditions: 



geneous fc-core percolation. The most interesting critical 
phenomena involve mixtures which separately give origin 
to phase transitions of different order. In this set of cases, 
we observe two main phase diagram topologies: a critical 
line intercepted by a first order line ending in a critical 
point, or a critical line matching the first order line onto 
a tricritical point (TCP). A careful analysis of equations 
Q and ([5]) for different types of binary mixtures shows 
that the occurrence of a TCP scenario is strictly related 
to the relative order of magnitude of (k a + l)-shells and 
the fc(,-corona (when k a < kf,) in approaching the hybrid 
transition. Broadly speaking, this condition quantifies the 
closeness of the parameters k a and kt in order for a TCP 
to occur. This criterion may be potentially important 
in network engineering (for example in managing large 
scale infrastructures), as the presence of a TCP implies a 
smooth change in the way the network collapses, and the 
possibility to control the order of the transition. 



h lk (0,Z) < f lk (Z) 



(52) 



Substituting the second equation into the first one yields 



+ OG 



q(q-l)P(q) 



q=2 W 

+(i---)£ ? Sri:( , 7 1 ) !Z '< 1 - z » , ^' so 

q=k W/ l = k V / 

(53) 

which is never true, because both sums on the left hand 
side are strictly positive for any Z > 0. 



It is interesting to note that the expansion (47) has 
terms which are reminiscent of a /c-shell decomposition 



(as in Eq. (|42j)). The violation of the TCP condition ([53} 
appears to be related to a mis-match between 2-shells of 
1-nodes and the leading term in Z of fc-nodes. 

Both in the cases k = (2, k) and k = (1, k), then, we 
can say that a necessary condition for the occurrence 
of a TCP is the vicinity of the two values of k. In the 
case k = (1, k), a TCP never occurs, whereas in the case 
k = (2, k) it is only present for k = 3, where the 3-shell 
of type 2 vertices has the same order of magnitude as the 
corona of type 3 vertices. 



This behavior is also reminiscent of the scenario ob- 
served in binary mixtures of hard spheres, where a dis- 
continuous glass-glass transition, ending into a critical 
point, weakens up to disappearing when the sizes of the 
two particle types become similar [33] . The analogy is 
even more striking, as recently proposed facilitated spin 
models, which belong to the same universality class of 
this HKC model [53], show with increasing evidence the 
correct reproduction of high-order singularities of the 
mode-coupling theory [32] and a strong link with kineti- 
cally constrained models [33]. In other words, high and 
low ks in fc-core percolation may be associated with large 
and small particles, respectively, as a large particle re- 
quires on average more neighboring particles to become 
caged and vice- versa for small particles. 



In summary, heterogeneous k-cove percolation provides 
an exactly solvable model characterized by a wealth of 
interesting critical phenomena. The model has a great 
potential in investigating the fundamental processes oc- 
curring in the proximity of critical points and therefore in 
giving insights into quite different interesting applications. 



V. CONCLUSIONS 

In this paper, we have explored and classified the critical 
phenomena which characterize binary mixtures in hetero- 
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